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Abstract: A general class of cosmological models driven by a nonlocal scalar field inspired 
by string field theories is studied. In particular cases the scalar field is a string dilaton or a 
string tachyon. A distinguished feature of these models is a crossing of the phantom divide. 
We reveal the nature of this phenomena showing that it is caused by an equivalence of the 
initial nonlocal model to a model with an infinite number of local fields some of which 
are ghosts. Deformations of the model that admit exact solutions are constructed. These 
deformations contain locking potentials that stabilize solutions. Bouncing and accelerating 
solutions are presented. 
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1. Introduction 



Field theories which violate the null energy condition (NEC) @, § are of interest for the 
solution of the cosmological singularity problem and for the construction of cosmological 
dark energy models with the state parameter w < — 1. 

One of the first attempts to apply string theory to cosmology Q was related to the 
problem of the cosmological singularity Q . A possible way to avoid cosmological singulari- 
ties consists of dealing with nonsingular bouncing cosmological solutions. In these scenarios 
the Universe contracts before the bounce . Such models have strong coupling and higher- 
order string corrections are inevitable. It is important to construct nonsingular bouncing 
cosmological solutions in order to make a concrete prediction of bouncing cosmology. 

Present cosmological observations || do not exclude an evolving dark energy (DE) 
state parameter w, whose current value is less than —1, that means the violation of the 
NEC (see |], [?J for a review of DE problems and || for a search for a super-acceleration 
phase of the Universe). 

A simple possibility to violate NEC is just to deal with a phantom field. The phantom 
field is unstable. There are general arguments that coupled scalar-gravity models violating 
the NEC are unstable ([0-0] and refs. therein). At the same time a phantom model 
could be an approximation to a nonlocal model that has no problems with instability [|l4j . 
A simple example of such a model is a model with a scalar nonlocal action (e~^ 9 (j)) 2 . In the 
second order derivative approximation e~^ 3 ~ 1 — D g this model is equivalent to a phantom 
one but does not have problems with instability. This type of models does appear in String 
Field Theory (SFT)(see [15| for a review) and in the p-adic string models [16]. The model 
with the particular kinetic term mentioned above is in fact realized in the p-adic string 
near a perturbative vacuum and is expected to be realized in the Vacuum String Field 
Theory (VSFT) Q. 

The purpose of this paper is the study of this type of models. As a model we con- 
sider a SFT inspired nonlocal dilaton action. Distinguished features of the model are the 
invariance of the action under the shift of the dilaton field to a constant as well as a pres- 
ence of infinite number of higher derivatives terms. A more general family of nonlocal 
models loosing the invariance of the nonlocal dilaton is also considered. For special val- 
ues of the parameters the models describe linear approximations to the cubic bosonic or 
nonBPS fermionic SFT nonlocal tachyon models, or p-adic string models [14], [ 18]- [ 24] . 
The NonBPS fermionic string field tachyon nonlocal model has been considered as a candi- 
date for the dark energy Several string-inspired and braneworld dark energy models 
have been recently proposed (see for example ^5|-[^7]] an d refs. therein). About a study 
of the tachyon dynamics with the Born-Infeld action see |2^, |2£| [3C| ]. 

We discuss a possibility to stabilize the model that violates the NEC in the flat space- 
time at the cost of adding extra interaction terms in the Friedmann background. One of 
the lessons from a study nonlocal dynamics in the flat case is a sensitivity of the stability 



problem to the form of the interaction term [31, 32, 33, gj, 35, 18, pq] . We use the 



Weierstrass product to present the nonlocal field in terms of an infinite number of local 
fields 1 22]. Some of these local fields are ghosts, which violate the NEC and are unstable. 
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The model is linear and admits exact solutions in the flat space-time. In non-flat case we 
get the same exact solutions after a deformation of the model. We used a similar approach 



to construct effective SFT inspired phantom models [37, 18, 35]. 

Another recently proposed model which violates the NEC and has higher derivatives 
is the ghost-condensation model [pSQ . Vector-scalar and tensor-scalar models that violate 
NEC and are stable in some region have been proposed in [4C, 41], respectively. 

The paper is organized as follows. In Section 2 we describe our strategy to the study of 
stringy inspired models. In Section 3 we present general solutions of the models in the flat 
case. Then we use some approximation to study these dynamics in the Friedmann metric 
and discuss cosmological properties of the constructed solutions. 



2. Set up 

In this paper we consider a model of gravity coupling with a nonlocal scalar field which 
induced by strings field theory 



A i — 
a xy—g 



</>F(-D 9 /M 2 )«/> 




(2.1) 



M p is a mass Planck, M s is a characteristic 



where g is the metric, EL 



9 ~ V=9"» _ 
string scale related with the string tension a', M s = I/t/c/, is a dimensionless scalar 

field (tachyon or dilaton) , g^ is a dimensionless four dimensional effective coupling constant 

related with the ten dimensional string coupling constant go and the compactification scale. 



A 



Mi 
in 



A! is an effective four dimensional cosmological constant. 



The form of the function F is inspired by a nonlocal action appeared in string field 
theories. In particular cases 



F(z) 



-fz + 1 



ce 



-2z 



(2.2) 



£ is a real parameter and c is a positive constant. Using dimensional space-time variables 
and after a rescaling we can rewrite (^^) for F given by (|2.2| ) as follows 




A' 



(2.3) 



where 



with the Plank mass [21 



and nip 



g^Mp/Mg. Generally speaking the string scale does not coincide 



] . This gives a possibility to get a realistic value of A |21[] . 
The form of the term (e^ 9 (j)) 2 is analogous to the form of the interaction in the action 



for the string field tachyon in non-flat background [14], which is a generalization of the 

. At some particular 



15, 32, 



SFT tachyon interaction term in a flat background [42, 

values of £ 2 and c this action appears in a linear approximation to SFT actions |44 
[48] and in a non-flat background has been considered in [14, 18, |l^, 20, 24]. The case 
of the open Cubic Superstring Field Theory (CSSFT) |l4], [2^] tachyon corresponds to 
£ 2 = - 1/ (a In (^g)) ~ 0.9556 and c = 3. We consider in detail action (|]| at c = 1, 
which is invariant under translation 6 4> + const. 
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We take the metric in the form 

ds 2 = -dt 2 + a 2 {t) [dx\ + dx\ + dxf) (2.4) 
and get the following equation of motion for the space homogeneous scalar field (p: 

F(-V)<p = 0, (2.5) 

where 



V = - d 2 - 3H(t)d t , H=- and a = d t a. (2.6) 

a 



The Friedmann equations have the following form 



3H 2 — — ^ £, 
3H 2 + 2H = 5 V, 



(2.7) 



where the energy and the pressure are obtained from the action fl2.1|) using standard formula 
V = , V = dmg{£ : V,V : V}. (2.8) 



For the case of F given by fl2.2|) the energy and the pressure have additional nonlocal terms 
£ n il and £ rd2 |3|, ||, |(J 

£ = £fc + £p + £nl2 + £nll + A', 
V = £k — £p + £ n l2 — £nll — A . 

Nonlocal term £ n n plays a role of an extra potential term and £ n ii2 a role of an extra 
kinetic term. The explicit form of the terms in the R.H.S. of ( |2.9| ) is 

£ k = ^(d<P) 2 , £ p = -l{f-c{e v <j>) 2 ), 

Snll = cf (e^^) ( - Ve^<f>) dp, (2.10) 

£nl2 = -c f (de^<p) (de^J>) dp. 



Our strategy is the following: 

First, we study the dynamics of the model ( |2.3| ) in the flat case. 

• We show that solutions have the form of plane waves. There are special values 
of parameters for which plane waves should be multiplied on linear functions. We 
calculate the energy and pressure on the corresponding solutions. 

• We present dynamics of the nonlocal model in terms of an infinite number of local 
fields |2^] and show that the energy and pressure densities of the nonlocal model are 
reproduced by the energy and pressure densities of the corresponding local models. 
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For this purpose we use the Weierstrass product representation for the function F in 



F(z) 



=./(*) 



n 



07, 



(2.11) 



where are complex numbers, and represent the fiat analog of ( |2.lD as 

Sfiat = \J d A x4>F{—U)<f) ~ ~ £ [e n ^„e^- D ) (□ + a* + c . c .] , ( 2 .12) 

where □ is the d'Alembertian in the flat space-time. 

• We consider approximated models obtained by a truncation of number of local fields. 

Then we consider the Friedmann Universe. There are two ways to study dynamics in 
the Friedmann metric: 

• One can use the found expressions for the energy and pressure in the flat case, E{(j)Q) 
and P(4>o), to calculate the corresponding Hubble parameter Hq, then using this 
Hubble parameter calculate a perturbation of the flat solution of equation of motion 
and so on: 



+ MHo) + .... H = Ho( fo) + Hiifo, fa) + ... 



(2.13) 



• One can search for deformations of the model that admit the same exact solutions as 
in the flat case and try to argue that the deformed models describe the initial model 
with a good accuracy. 

Both ways permit to find the first approximations to the models (|2.5|) , ( |2.7D . The first 
way have been used in [20]. In this paper we will follow the second way. 

To this goal we use a representation of nonlocal dynamics given by action ( |2.lD in 
terms of local fields 

,2 



J <fx^ I + [^ne f{ - D ^(n g + a 2 n )i> n + A + c.c. 



(2.14) 



We perform a deformation of this model by several steps. First, we consider an approxi- 
mation to the model ( 2.14j ) in the form 



/ d'xV^g {^R + £ [^Vw^D, + <4)V>n + A' + cc. 



(2.15) 



Second, we restrict a number of local fields and, third, we add potentials of the order 1/rrip 
in which A' is also included: 



/ d'x^-g {^R + £ [|^ne^(n fl + al)^ n + c.c 



V(tpl,...,tp n ) 



(2.16) 



such that solutions of the field equations in the non-flat case are the same as the flat case. 

Finally, we find the corresponding scale factor a(t) and study cosmological properties 
of approximated solutions to our model. 
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3. Flat Dynamics 
3.1 General Solutions 

3.1.1 Roots of the Characteristic Equation 

In the flat case the action fl2.1|) has the following form: 



Sflat = \j d 4 x0F(-D)«/». (3.1) 

Equation of motion on the space-homogeneous configurations ( |2.5| ) is reduced to the fol- 
lowing linear equation: 

F{d 2 )(/) = 0. (3.2) 

A plane wave 

<j> = e at (3.3) 



is a solution of (|3.2|) if a is a root of the characteristic equation 

F{a 2 ) = 0. (3.4) 



For a case of F given by (2.2) equation (|3.2| ) has the following form 



-fd 2 (j) + (j) - ce- 292 (/) = 0. (3.5) 
This equation has an infinite number of derivatives and can be treated as a pseudodiffer- 



ential as well as an integral equation [35|. 



The corresponding characteristic equation: 

F{a 2 ) = -fa 2 + 1 - ce~ 2a2 = (3.6) 

has the following solutions 



1 / 2ce~ 2 /£ 2 » 
a n = ±-J4 + 2eW n ( _j, n = 0,±l,±2,... (3.7) 

where W n is the n-s branch of the Lambert function satisfying a relation W(z)e w ^ = z. 



The Lambert function is a multivalued function, so eq. (3J3) has an infinite number of 
roots. 

Parameters £ and c are real, therefore if a n is a root of ( |3.6|) , then the adjoined number 
a* is a root as well. Note that if a n is a root of (|3.6|) , then — a n is a root too. In other 
words, equation ( |3.6| ) has quadruples of complex roots 

a n ,±± = ±Re(a n ) ± ilm(a n ). (3.8) 

If a 2 = a.Q is a multiple root, then at this point F{afy = and -F'(qq) = 0. These 
equations give that 

2 1 1 . . 

«0 = 71 - 2' ( 3 - 9 ) 
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hence ckq is a real number and all multiple roots of F{a 2 ) ) = are either real or pure 
imaginary. The multiple roots exist if and only if 

Real roots for any £ and c, except £ 2 = and c = oo, are no more then double 
degenerated, because F"{otfy / 0. 

Summing up we note that according as the values of parameters c and £ 2 there exist 
the following types of the general real solution of (|3.5|): 



If c 7^ l~e 2// £ 2 and c 7^ 1 then the general real solution is 



n n 

where i? n and C n are arbitrary real and complex numbers respectively. 

If c = ^e 2 ^ 2 > 1, then to get the general real solution one has to add to ( |3.11| ) 



>o-ff./< ' + R->I<- '. if e 2 < 2, (3.12) 

>„ - C ',/,'""' + (:"',/<" '■ "n-/ A /^-^ if e 2 > 2. (3.13) 



If c = 1 then to get the general real solution one has to add to ( |3.11 ) 



0o = Cit + q> s if f^2, (3.14) 
0o = C 3 i 3 + C 2 t 2 + dt + C , if £ 2 = 2. (3.15) 

3.1.2 Real Roots of the Characteristic Equation 

For some values of the parameters £ and c eq. ( |3.6| ) has real roots. To mark out real values 
of a we will denote real a as m: m = a. 

To determine values of the parameters at which eq. ( |3.6| ) has real roots we rewrite this 
equation in the following form: 

e 2m 2 _ c 

£ 2 = g(m 2 ,c), where g(m 2 ,c) = — - 2 . (3.16) 



ra'e 



The dependence of g(m, c) on m for different c is presented in Fig. ([!]). This function has 



a maximum at 



2 1 1~ / e 



-1 



m 2 max = --- gWLi ( — ) , (3.17) 



provided c is such that W_i < — 1; i* 1 the other words < c < 1. 

There are three different cases (see Fig. [l]). 

• If c < 1, then eq. (|3.6|) has two simple real roots: m = imi for any values £. 
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If c = 1, then eq. ( |3,6| ) has a zero root. Nonzero real roots exist if and only if £ 2 > 2. 
If c > 1, then eq. ( |3.6| ) has 

— no real roots for £ 2 > ^ ax , where 

1 — ce -2 " 1 " 1 ^ 2 



ra L W_i(-e 7c) 



(3.18) 



— two real double roots m = ±m max for £ 2 = £^az 

— four real single roots for £ 2 < Cmax ■ ^ n this case we have the following restriction 
on real roots: m 2 > ^ In c. 

3.1.3 Pure Imaginary Roots of the Characteristic Equation 

For some values of the parameters £ and c eq. ( |3.6| ) has a pair of pure imaginary roots. Let 
us introduce a new real variable /i = ia. From ( |3.5| ) we obtain 

(£V + l)e- V = c. (3.19) 

This equation is equivalent to 

c _ e - V 

£ 2 = <?0 2 >c), where g(p 2 ,c) = — — 5-5-. (3.20) 

/i z e z ^ 

The dependence of g(fJ, 2 , c) on // for different c is presented in Fig. ||. 
For different c we have: 

• For c < 1 there are two real simple roots \jl = ±//i, 

• For c = 1 nonzero real roots exist only if £ 2 > 2, 

• For c > 1 real roots exist if and only if 



£ £,min — / TV (3.21) 
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0,5 1 1,5 A* 2 0,5 1 1,5 ^ 2 



Figure 2: The dependence of function g([i 2 , c), which is equal to £ 2 , on /i at c = 1/2 (right), c = 1 
(center) and c = 2 (left). 



If £ = Cmin, then there exist two double real roots: fi = ±fi m i n , where 



Umin = 2 + 2W (~~~ ) ' 



(3.22) 



At £ > £ TO mj eq. ( 3.19| ) has four real simple roots. 
3.1.4 Roots at the SFT inspired values on £ 2 and c 

Let us consider a special values of £ 2 and c, which have been obtain in the SFT inspired 



cosmological model [14, 20 1. From the action for the tachyon in the CSSFT |45|] the 
following equation has been obtained 



(-e a 2 + 1) = 3e" 52 / 4 , 



where 



f 2 



0.9556. 



(3.23) 
(3.24) 



Substituting a = 2y/2a, we obtain eq. ( |3.6| ) with £ 2 = 8£q and c = 3. From ( |3.10| ) it is 
follows that all roots are simple. We obtain that £^ in > £ 2 > ^ a3; , so there exist neither 
real roots nor pure imaginary roots. 

3.2 Energy Density and Pressure 
3.2.1 General Formula 



Equation (3.5) has the conserved energy (compare with |4l], |5(], jig]), which is defined by 
the formula that is a flat analog of (|2.9| ). The energy density is as follows: 



where 



Ek 



E — Ek + E p + E n n + E n i2, 



■my, 



e p = Af + |* 2 , 



(3.25) 
(3.26) 
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E nll = c J' (e^*) (d\e» 92 <S>)) dp, E nl2 = -c j' (<9( e -^ 2 <I>)) (d(e^ 2 <I>)) dp. 

(3.27) 

For the pressure 

P = E k — E p — E n n + E n i2, (3.28) 



we have the following explicit form 



P = ^(S0) 2 + 



2 v T/ 2 



^2_£$a_ c f'^-pd*^ (sP^^)) _ (a(e-^ 2 $)) (d{e? 92 <$>))} dp. 

(3.29) 

Let us calculate the energy density and pressure for the following solution 



A? 



°ne ant , (3.30) 



n=l 



where N is a natural number, C n are some constant and a n are solutions to eq. (|3.q ). 
For N = 1 and 

= Ce a * (3.31) 

we obtain 

£(Ce Q *) = 0, (3.32) 

P(Ce at ) = C 2 p a e 2at . (3.33) 

Hereafter we denote the energy density and pressure of function (p(t) as the functionals 
E((f)) and P(4>), respectively, and use the following notation 

Pa = « 2 (e 2 - 2 + 2£V) . (3.34) 

For JV = 2 and 

= C ie ait + C 2 e a2t , (3.35) 
where a\ and «2 are different roots of (|3.6| ) we have (see Appendix A for details) 

E (Cie Ql * + C 2 e" 2 *) = -2CxC 2 p ai , at a 2 = -ai, (3.36) 

and 

E (C ie ait + C 2 e a2t ) = 0, at a 2 ^-a x . (3.37) 
The pressure P(0) for solution ( |3.35| ) is 

P (C ie - at + C 2 e at ) = (C 2 e~ 2at + C\e 2at ) p a . (3.38) 

In the general case we have 

/ N \ N N 

E\Y^C n e a A = -2^ £ C„C fe Pa n <5a„,- Qfc , (3.39) 

\n=l / n=l fc=n+l 
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where 

j 1, a n = -a k , 
S an ,-a k = < (3.40) 
^0, a n ^ -a k . 

Note that all summands in ( p. 39 ) are integrals of motion, therefore, we explicitly show 
that E (^Yln=i C n e ant ^j is an integral of motion. From formula ( |3.39| ) we see that the 
energy density is a sum of the crossing terms. At the same time the pressure is a sum of 
"individual" pressures and has no crossing term. In the case of an arbitrary finite number 
of summands the pressure is as follows: 

/ N \ N N 

P £ C n e^ \=Y, C l P = E c2 nP^e 2ant - ( 3 - 41 ) 

\n=l / n=l n=l 

If the parameters and c are such that the characteristic equation ([O]) have double 
roots, then eq. ( |3.5| ) has the following solution 

Mt) = B ie a °H + de aot + B 2 e- aot t + C 2 e~ aot , (3.42) 

where B\, B 2 , C\ and C 2 are constants, Qo 7^ is defined by (|3.9|), Using formulas ( 3.25| ) 
and ( p.28| ) and substituting 



«o = 2g » ( 3 - 43 ) 



we obtain 



E(M = ~ ^ 2 ^ (3Ci5 2V /4-2^ - 3C 2 Siv/4-2^ + 8S X S 2 (2£ 2 - l)) . (3.44) 
The pressure is as follows 

P(0o) = - V 3 ^ 2 ^ ( 5 2 e — (B 2 (Sf - 3^4-2^ - 4j - 3^^4-2^) + 

+ ^e^? 2 ^3CiV4-2^ + Bi (8£ 2 + 3*^4 - 2£ 2 £ " 4 ))) • 

(3.45) 



3.2.2 Energy Density and Pressure for real a 

As we have seen in Sect. 3.1.3 for some values of parameters £ and c eq. ( |3.6| ) has real 
roots. We denote as p m the values of p Q for real a = m, 

Pm = m 2 (e 2 - 2 + 2£ 2 m 2 ) , (3.46) 

where £ 2 is given by ( |3.16| ) . For different values of c the function p m is presented in Fig. |3|. 

If and only if c > 1, then there exists the interval of < m 2 < m max on which 
p m < 0. Some part of this interval is not physical because on this part g(m 2 ,c) < 0. The 
straightforward calculations show that 

_ dg(m?,c) 

Pm- g m2 ■ (3.40 
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Figure 3: The dependence of p m on m at c — 1/2 (right), c = 1 (center) and c = 2 (left). 
Therefore, at the point 



2 2 



j2 

^max 



(3.48) 

we have two positive 



we obtain p m (m mo3; ) = 0. We conclude that for c > 1 and £ 2 < ^ 
roots of (|3~lj|): mi and m-i > mi, with p mi < and p m2 > 0. 

The energy density and the pressure for solutions with real a one can calculate using 
formulas ( |3.39| ) and ( [3.41 ) and results are presented in Table 1. 





E 


P 







±2mt 

Pm c 


sinh(mi) 




p m cosh(2mt)/2 


cosh(mi) 




p m cosh(2mt)/2 



We see from Table 1 that odd solutions 
are physically meaningful, E > 0, if p m is 
positive. Fig.|| shows that odd solutions are 
physical for c < 1 and any £ 2 and for c > 1 
only for m 2 > m 2 max . The pressure corre- 
sponding to this solution is always positive. 

Even solutions are physically meaningful 
if Pm is negative. Therefore, even solutions 

are physical only for c > 1 and m 2 < m max . The pressure corresponding to this solution is 
always negative. The equation of the state parameter for this solution is 

w = — cosh(2mi) < — 1. 

3.2.3 Energy Density and Pressure for pure imaginary a 



Table 1: Solutions, densities of energy and 
pressures for real a 



As we have seen in Sect. 3. 1.4 for some values of parameters £ and c eq. ([O 
of pure imaginary roots. These solutions correspond to 

4> = C sin( / u(t — to)). 

On the solutions (3.50) the energy and pressure are given by 



C 2 C 2 
E(Csm(fi(t - t ))) = —p a = -yTT/i, 

C 2 

P{Csm(p(t-t ))) = -— 7r M cos(2/i(t-t )), 



where 



7T„ 



= fi (2+ 2ev - e) 



(3.49) 

has only pair 
(3.50) 

(3.51) 
(3.52) 

(3.53) 
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where £ 2 is given by ( 3.20| ). tt^ as a function of /i for different values of c is presented in 
Fig (1). 




0,2 0,4 0,6 




0,2 0,4 0,6 0,8 // 1 




Figure 4: The dependence of ^ on /i at c = 1/2 (right), c = 1 (center) and c = 2 (left). 

Note that ir^ for c ^ 1 is positive, however for c > 1 is positive only for fi 2 ^ Mmin 
| — p . The equation of the state parameter for this solution is 



w = — cos(2fj,(t — to)), = 

3.2.4 Energy and Pressure in the case c = 1 

The energy density and pressure for solutions ( [3.14 ) 



\w\ 1. 



(3.54) 



0i (t) = dt + C , 
where Cq and C\ are arbitrary constants, are as follows 

mi) = (f - 1) Ci 2 , p(0i) = (|-i)(f 



(3.55) 



(3.56) 



and the state parameter w = P/E 1 = 1. 

The straightforward calculations show that pressure and energy density for more gen- 
eral solutions 



A' 



J2Cne ant + fa 



(3.57) 



n=l 



arc 



and 



E (^2 Cne ant + fa^j=E Cne a "^j + E (fa) (3.58) 



N 



\n=l 



+ P(fa) 



(3.59) 



Let, for example, 



fat) = cosh (mi) + Cit + Cq. 
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The corresponding energy density and pressure are: 

E = -\v m + (y - l) Cl P = ^ m cosh(2mi) + - 1 ) ( f . 



In the case £ 2 = 2 the root a = is a double root of eq. (p.6|), so eq. ( |3.5| ) has solu- 
tions (l3~ll ): 

02 = C 3 i 3 + C 2 t 2 + Cit + C7o. 

We obtain: 



£(^ 2 ) = 4 (3dC 3 - C 2 - 6C|) , P(0 2 ) = 72C|t 2 + 48C 2 C 3 t + 4 (C 2 + 3dC 3 - 6C|) . 

3.2.5 Energy and pressure for complex a = r + iv 

For a decreasing solution 

cf)(t) = e~ rt cos(i/i) (3.60) 

we have 

E 1 (e~ r * cos(z^)) = 0, (3.61) 

— 2rt 

P{e- rt cos(vt)) = (p ae 2i^ +pa , e -2^) (362) 

Using general formulas (3.39) and (|3.41| ) it is easy to find the energy and pressure for 
even and odd real solutions with complex a. For example the energy and pressure for the 
even solution <p = cosh(rt) cos(i4) are as follows: 

£7(cosh(ri) cos(i/t)) = -- (p a + Pa*) , (3.63) 

o 

P(cosh(rf) cosfi/t)) = V — (e~ 2at + e 2at ) + ^ ( e" 2 "*' + e 2a * 1 ) . (3.64) 
16 16 V / 

The equation of state parameter on the even solution is 

Pa + Pa /n C r\ 

w = 1 7 i \ I ; *,y {3.65) 

p a cosn(ar) + p* a cosn(a*tj 

3.3 Local Field Representation 

3.3.1 Weierstrass product and Mode Decomposition for the Action 



As in [22] we can present F{— □) in the action (2.1) as the Ostrogradski Representation. 
To this purpose let us construct the Weierstrass product for the function F(z) of a complex 
variable z. Let us recall that a complex function R(z) such that its logarithmic derivative 
R'(z)/R(z) is a meromorphic function regular in the point z = 0, has simple poles and 
satisfies \R'{z)/R(z)\ < C, z £ T n , n = 1,2,..., can be presented as 

R(z) = R(0)e^m z ]J I 1 - — J e z ' Zk . (3.66) 

T n , n = 1,2, ... is a set of special closed contours T n such that the point z = is in all 
T n , T n is in T n+ i, and S n /d n < C, where S n is a length of the contour T n , and d n is its 
distance from zero [51|. 
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In the case of a more week requirement \R'(z)/R(z)\ < M\z\ p z £ T n , n = 1,2, ... we 
have 

R(z) = II " |) ^ QkM =Ej(^)'' ( 3 - 67 ) 

where /(z) is an entire function. 

In the case ol R = F given by fl2.2| ) the Weierstrass product can be written in the form 

F(a 2 ) = e f(a2) Y[(a 2 -al). (3.68) 

n 

The function f(z) in our case is 

f(z) = A + (3z, (3.69) 

where constants A and /3 are determined by £ and c. It will be shown that the equations 
of motion do not depend on values of A and (5. 

It is convenient to pick out real roots in (3.68) and combine the complex conjugated 
roots: 

F{* 2 ) = e A+ ^ 2 H(a 2 - ml) J](« 2 - a 2 n )(a 2 - < 2 ), (3.70) 

where rrik denote real roots. In Subsection 3.1 we have found the cases when real roots do 
exist. 

In the case of simple roots the Lagrangian up to a total derivative can be presented as 
a sum of an infinite number of fields [ 52 ] — 1 55 ] 

C = ^F(d 2 )0 ~ i £ \e n ^ d2 \-d 2 + a 2 n )^ n + c.c] , (3.71) 

where ~ means equivalence up to a total derivative, e n are constants. It is the Ostrogradski 
representation. Note that for complex roots ip n are complex. 

3.3.2 Mode Decomposition for Energy Density and Pressure 

It is instructive to present the formula for energy and pressure obtained in section 3.2. in 
terms of tp fields. All considerations below take place for nondegenerate roots. 

According to a general procedure of construction of the energy and pressure we write 
a generalization of ( |3.71| ) to a non-flat case 



C g = Cg^ n ), C g {^ n ) = ^^Vne /( ~ D9) (□<, + a 2 ) </> n , (3.72) 



and find 



= J>n, E n = ^ (>„ 2 - a 2 j£) (3.73) 

n 

P i> =Y J P n, Pn = ^ Un + O^ 2 ) . (3.74) 



The E.O.M. for ip n is 



{3 2 - a 2 n )i{j n = (3.75) 



- 15 - 



and its solutions are 

^ n = A n e a - t + B n e- a - t . (3.76) 

For solutions ( [3. 76 ) we obtain 

= 2 a 2 n e n A n B n e^ , (3.77) 

n 

P i) = ^2 tnc? n {A 2 n e 2a ^ + Ble' 2 ^) e** . (3.78) 

n 

On the other hand according to ( |3.39| ) and ( |3.41| ) we have 

E = _ 2 ^2 A n B n a 2 n (f - 2 + 2fa 2 n ) , (3.79) 

n 

P = E {^~ 2ant + By ant ) < (e 2 - 2 + 2iW n ) . (3.80) 

n 

Comparing ( |3.77| ),( |3~7"8 ) and ( p. 79P , ( 3~80| ) and using equation fl3.6j ) we obtain that 

E = Ej/j, P = P^ (3.81) 

if and only if 

e n = -(2ce- 2a *+f)e- f3a *, (3.82) 

that is in accordance with general formula for e n |34|, ^] . Note that we consider only simple 
roots a n . 

3.4 Finite Order Derivative Approximation 
3.4.1 Two types of approximations 

There are two different types of finite order derivative approximations: 

• a direct finite order derivative approximation 

k 

L 2k (a 2 ) = Y,a n a 2n ; (3.83) 

n 

• an approximation by a finite number of terms in the Weierstrass product 

L$P (a 2 ) = f(a 2 )X{(a 2 -a 2 n ). (3.84) 

n 

We label roots so that \a n \ ^ |a n+ i|. 

Locations of roots of the characteristic equation in the complex a-plane are presented in 
Figure |5[ One can see that structures of roots location for different values of £ look similar. 
Let us consider the most simple case £ = 0, c = 1, that corresponds to 

L(a 2 ) = L ,i(a 2 ) = 1 -e" 2 " 2 . (3.85) 
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In this case all roots can be written explicitly 

„4 



1 



-2a* 



2a 2 e 



n 

3=1 



i + 



or 



(3.86) 



Note that all conclusions can be generalized on the 
case of arbitrary £ and c such that all roots a n are 
simple. 

3.4.2 Direct Finite Order Approximation 

In the second order derivative approximation we sho- 
uld keep in (3.5) with c = 1, £ = only the second 
derivatives 





4- 








+ 

- + 

+ + + + 2 : 
+ 


+ 
+ 


+ + 
+ 


-4 


2 4 




2 ] 




+ 

- + 


+ 
+ 


+ + 










-P+ + 

-H+ + 


-4 




-t r +H- 



20 2 



0. 



(3.87) 



This equation has the following solutions 



Cit + C 2 . 



Figure 5: Roots for c = 1 and £ = 

(big crosses), £ = 2 (middle crosses) 
and £ = 15 (small crosses). 



(3.88) 



These solutions correspond to roots ot\ = 0. 

In the fourth order derivative approximation eq. (3.6) at c = 1 and £ = is as follows 



2af - 2a 4 







(3.89) 



and has two solutions: 

a 2 = 0, a 2 = 1. (3.90) 

We see that in the fourth order direct approximation the approximate equation ( |3.8£| ) 
has a root (a = 1) that is absent in the initial equation ( |3.6| ). 

The similar situation takes place for higher order approximations. The characteristic 
equation of the direct n-order approximation contains several artificial roots. The appear- 
ance of these roots is related to artificial roots of polynomial approximations of the function 
/ in the Weierstrass product ( 3.68[) . In Figure |6| we plot all roots of an n-order polynomial 
approximation of function Lo,i( Q;2 ) f° r n=20 and 40. We see that these polynomials have 
true roots as well as artificial roots that go to infinity when the order of the approximate 
polynomial increases. 

3.4.3 Finite Order Weierstrass Product Approximation 

We consider an approximation to -Lo,i(a 2 ) that keeps only a finite number of terms in the 
Weierstrass product 

k 



L^ w) ( a 2 ) = 2a 2 e~ a2 

^0,1,2+4^° > — z " e 



n 

3=1 



1 + 



a 



n 2j2 



(3.91) 



We call this approximation as (2k+l)-roots approximation. The corresponding Lagrangian 
is 



C 



(W) 

0,1,2+44 



1 



(3.92) 
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1 2» 



A. 



B. 



C. 



Figure 6: A. Roots of the 20-th order approximation Lo,i;io( a ) of the function Lo,i( a2 ) depicted 
by small diamonds. There are one root in the origin, 4 roots located on a "cross" and 15 roots on 
two bow-like curves. B. Roots of the 40-th order approximations, Lo,i;2o(a0, and roots of -Lo,i;io(<*) 
depicted by small boxes and small diamonds, respectively. £o,i,2o(o0 has one root in the origin, 8 
roots located on a "cross" and 4 of them coincide with true roots of -Lo,i;io( a )- Artificial roots of 
A),i,2o(a 2 ) are located on bigger "bows" as compare with artificial roots of Lo,i,io(a 2 ). C. Exact 
roots of Lo t i(a) are located on a cross and depicted by small crosses. The roots of the L ^. w {a) 
and Lo,i;2o(o ; )) are depicted by small diamonds and small boxes, respectively, as in A and B. 



Let us consider one root approximation 



45 - 



E.O.M. is 



2-d 2 



d z e 







(3.93) 



(3.94) 



and it has an unique solution (the same as E.O.M. of 2-nd order derivative approximation) 



<f> = At + C. 

The next approximation is 3 roots approximation 



£SJ = i*fl»«-*(i + £)* 



E.O.M. is 



d 4 



d z 1 + -z )e-° 4> = & d z 1 + -= U = 



d 4 



7T Z 



and it has the following solutions 

<f>e(t) = A 1 t + S 1 + A 2 e^ t sm(^t + d2) + A 3 e-^ t S m(^t + S 3 ), 

where A^ and 5^ are arbitrary constants. 

The next approximation contains 4 extra roots 



4vr 2 



(3.95) 



(3.96) 



(3.97) 



(3.98) 



(3.99) 
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E.O.M. is 

52 ( 1 + S)( 1 + S) e ~^ = ° (3 - 100) 
and it has the following solutions 

lo (t) = <fo(t) + A 4 e^ 1 sm(^t + 64) + A 5 e-^ 1 sm(^t + S 5 ). (3.101) 

It is obvious that solutions of the Weierstrass approximate equations reproduce a finite 
number of modes of the full equation. 

If we restrict ourself to decreasing solutions 

^io,decr(*) = *i + MeT^ 1 sin t + S^j + A^^ 1 sm(^t + 6 5 ) (3.102) 

we see that the last term can be ignored as compared with term first two terms. 

We can see that the Weierstrass product approximation is more preferable than direct 
approximation, because there is no problem with extra roots. In the next section we will 
use the Weierstrass product approximation to construct local cosmological models. 



4. Non-flat Dynamics 



4.1 Modified Action 



The goal of this section is to consider the nonlocal model (2.1) in the Friedmann Universe. 
To consider the dynamics in such a system we need to solve nonlinear Friedmann equations 
fl2.7| ), which represent hopelessly complicated problem. From (|2.7| ) we obtain the following 
nonlinear integral equation in H(t): 



H 



2mj; 



(V + £) 



-2 (t W) 2 

nip \ 2 



de^ v 4>) (de^-PV^ dp J , (4.1) 
3H{t)dt. There is no method to solve eq. (|4.1|), even if the function <p{t) 



where T> 1 
is given. 

As a first step to consider such a problem we would like to construct some effective 
exactly solvable models that can be consider as an approximation to the nonlocal model. To 
do this we use the finite order approximations, constructed by the Ostrogradski method. In 
other words we choose a special solution of eq. (3.5) and find the corresponding Ostrogradski 
approximation in the flat space-time. After we deform the obtained approximate model to 
the case of the Friedmann Universe, assuming that exact solutions in the Friedmann metric 
are coincide with exact solutions in the flat space-time. Note that the similar assumption 



has been used in the papers [37, 18, 38 1 to construct effective local models with exact 
solutions. 



Our starting point is the Lagrangian ( 3.71 ). The corresponding action in the non-flat 
space-time is as follows: 

n2 



So 



str. 



d A x, 



111 



-9 



[1^ 



(4.2) 
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If the fields <p n depend only on time and the metric is a spatially flat Friedmann metric, 
then W6 have the following equation for ip n 

e n (V + a 2 n ) i> n - e-^»V; n = & (2c e " 2a ' + £ 2 ) (V + a 2 n ) ^„ + = 0, (4.3) 

where VL is a derivative of V on Tp n . Note that form of V(tpi, ... ,ip n ) depends on choose 
of special solutions ipi, . . . , ifi n . The form of V(ipi, . . . , ip n ) is given below (Subsection 4.4). 
The energy and the pressure density in the Friedmann metric have the form 



£mod — + V, 
Pmod Pip 



(4.4) 
(4.5) 



where E^ and Pw, are given by formulas ( 3, 73] ) and ( |3.74 ) respectively. This means that 
the extra term V play a role of a potential term. 
The Friedmann equations of motion are 



3H 2 



3H 2 + 2H 



mt 



(E^ + V) , 
ml 



Therefore 



H 



P%p + Ej[,) . 



(4.6) 



(4.7) 



We choose such V that ip^ in the non-flat case are the same as in the flat case. Using 
fl3~77p and (jTF|)) we get 



k = ^2 J2 a l^ al (Ale 2 ^ + 2A n B n + B 2 e~ 

P n 

Using (3.76) we can rewrite (|4.8|) as follows 



■2a n t 



-2o? n + ^2\ .2 



(4.9) 



that 
H = 

where 



2ml ^ n 2ml ^\ 

P n P n 

Substituting values of e n (formula |3.82 ) and using formulas ( 3.39j ) and ( |3.41| ), we obtain 

i (e 



2m\ ^ 2 



(30) 2 - c£ (a(e"^ 2 $)) (d(e?**j) dp) = W) + (4-10) 



(4.11) 



Therefore, 



P V n n 



K e -2a n t + ^l e 2a n t ) p 1 + Hqj ( 4 . 12 ) 



2a 



2a r 
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where H$ is an integration constant and we assume the sum goes over the complex conju- 
gated roots. 

It is convenient to rewrite fl4.9| ) as follows 



H 



— Y 

2ml ^ 

P n 



Pa, 
al 



1 



2m 2 



(f - 2) 



n 



n 



2 12 



(4.13) 



Thus to obtain the crossing of cosmological constant barrier one should consider the 
case £ 2 < 2 and the field 4>(t), which consists of at least two modes. It is easy to see 
that H(t) has no singular point at finite time. For some values of parameters we obtain 
bouncing solutions, which satisfy the conditions H(0) = and H (0) > 0. 

In the following subsections we construct effective potentials for one-, two- and N- 
modes solutions. One-mode models can be consider as toy-models, whereas two-fields 
models are more realistic. Note that there are modern inflation models, for example [5(| 57], 
which include two scalar fields (see also [^] and references therein). In this paper we 
describe a procedure to construct effective models and analyse only the simplest properties 
of them. More detail analysis and compare with the nonlocal model dynamic is a subject 
of our future investigations. 



4.2 One Mode Solutions 

4.2.1 One Real Root: General case 



Let q is a real root of ( |3.6|) , without loss of generality we can assume a > 0. The corre- 
sponding scalar field is (see ( 3.76[) ) as follows 



ijj(t) = Ae at + Be~ at . (4.14) 
The function ijj(t) is a solution of the following first order differential equation: 

^ 2 = a 2 (^ 2 - AAB) . (4.15) 



Using superpotential method j59[ (see also |37|, |38|) we consider the Hubble parameter 
as a function of if): 



(4.16) 



From fl4~13| ) we obtain 



dW 



Pa 



2mj l a 2 



=F 



2m 2 ,a 



y/ip 2 - 4AB 



and 



W(ip) = ^^- (tViP 2 - 4AB - 2ABln (a ( V + v 7 ^ 2 - 4AB) ) ) + H . 



2ra^a \ 2 



The corresponding potential V(V>) is 



V(V0 = 3m 2 p W 



2w2 Vol >, 2 



2cr 



2w2 Vol , ,9 



{if 2 - 4AB) . 



(4.17) 



(4.18) 



(4.19) 
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We see that the potential depend on values of constants A and B, more exactly depend on 
value of the production AB, and does not depend on sign of W(ip). The potential V(ip) is 
not a polynomial, but in the limit of the flat space-time (m 2 — > oo) we obtain a quadratic 
polynomial. Note that in the case of the flat space-time one can eliminate constants A 
and B adding a constant to potential V(ip), so in this case the equations of motion do 
not depend on these constants. In the following subsections we analyse in detail a few 
particular cases. 

4.2.2 One Real Root: Decreasing Solution 

Let us consider a simplest particular case: 



From ( 4,12j ) we have 



Hit) 



tp = exp(— mt), 

_ Pm -2mt 

4nipm 



(4.20) 



(4.21) 



Therefore to ensure that ( 4.20| ) and ( 4.2 1| ) solve the Friedmann equations we have to add 
to the action the following potential 



v(V) 



lGnipm 



2 / 4 , 3p m 2 Pm , 2 , o 2 tt2 

■>Pm<' + ^ H °^ +3m p H . 



(4.22) 



Let us note that one gets the same potential for the unbounded solution ip\ = exp(mi). 
The decreasing solution ( 4.20[ ), ( |4.21 ) corresponds to the scale factor 



a(t) = ao exp 



8nipm 2 



exp(-2mt) + H t . 



(4.23) 



This solution has no singularities at finite t. It describes an increasing Universe only in the 
case p m < 0. As we have seen in Sect. 3.2.2] the negative pressure is possible if c > 1 and 



m < mi 



that takes place if £ 2 < S, max , where ^ m ax is given by ( 3.18 ). 



4.2.3 One Real Root: Odd Solution 

Let us consider a odd solution 

tp = sinh(mt). 

From (gTg ) we have 

H(t) = 



P m ■ urn j.\ Pm, , tt 

smh(2mt) — - — + Hq 



8m 2 m 



Am 2 



and the potential for Hq = is 



(4.24) 
(4.25) 

(4.26) 



The potential function of ip2 is given by 

ip 2 (l + tp 2 ) + 2V>\A + V ;2 arcsinh(V') + arcsinh(V') 2 . (4.27) 



vn>) = -?n + 3p2m 

2 lQwlm 2 
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The odd increasing solution ( |4.24| ), ( [4.25 ) corresponds to the scale factor 



a(t) = ciq exp 



Pr, 



16m 2 m 2 



(cosh(2mt) + 2m 2 t 2 ) + H t 



(4.28) 



The scale factor ( 4.28] ) has no singularity. It increases at large time if p m < 0. 

4.2.4 One Pure Imaginary Root 

Let us consider an odd solution 

tp = sin(/ii). 

From ( 3.51| ), ( |3.52| ) and ( [4.7] ) we get explicitly the Hubble parameter 

H{t) 



t + ~~~^~~t; sin(2/it) + Hq. 



Am 2 8/tm 

The potential for Hq = as a function of ip is given by 



(4.29) 



(4.30) 



v(V0 = ^ + 



3vr2 



2 " p ' 16m 2 {i 2 



ip 2 (l — ip 2 ) + 1 — ip 2 arcsin(V') + arcsin(-0)^ 



(4.31) 



Note, that this formula is valued only for ip 2 ^ 1. On this region the potential ( [4.31 ) is 
convex and it has an unique minimum. 

The odd periodic solution Q4.2S| ), ( [4.30 ) corresponds to the scale factor 



a(t) = clq exp 



16fj, 2 m 2 



cos(2/ii) + — ^t 2 + H t 



8m 2 



(4.32) 



and we have an expansion with an acceleration if 7T„ > 0. In this case at Hq = we obtain 
the bouncing solution. 



4.3 Solutions with the crossing of cosmological constant barrier 
4.3.1 Pair of Complex Roots 

For the case of a complex root a = r + iu we consider the following real solution 
ip4 = i>a + i>a* , where ip a = e~ at , 
From ( jO||) we have 



-a*t 



H(t) 



1 



4m 2 . 



Pa -2at , Pa* -2a* t , rr 

— e H — e +-no 

a a* 



4m 2 



Pa 1 2 , Pa* / ; * \2 , tj 

a a* 



(4.33) 



(4.34) 



The corresponding potential is as follows 
3 



V 



16m 2 



a 



if -2 + 2£ V) 4 iPi + a* 2 (e - 2 + 2£ V 2 ) ^ + H 2 



+ 2aa* (£ 2 - 2 + 2£ 2 a 2 ) - 2 + 2^V 2 ) ^ + 2tf + ^(C^ 



(4.35) 
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Prom QL34J) it follows that 



H 



2m 2 



-2at _,_ iP -2a*t 



-2rt 



2rra 2 



(Pa + Pa*) CO$(vt) +i{p a - p a * ) Sm(vt) 



It is easy to check that p a * = p* , so gi = p a + p a * and (72 = i(Pa — Pc 
numbers. Formula ( [4.36| ) can be rewrite in the following form: 



H 



2m 2 

where 7 is defined by the following relations: 



e Irt sm{vt + 7), 



sin(7) 



V<8 + ii ' 

As known the state parameter 



w 



■1 



cos(7) 



2 H 



72 



3H 2 ' 

so we obtain that w crossing the barrier w = — 1 infinite number of times. 
For example, in the case: £ 2 = and c = 1, when p a = —2a 2 and 



mi, £ — ' 



If we choose a = \pH + iy/ir, then p, 
1 



-2m, qi = 0, q2 = — 47r, 7 = 7r and 



mt 



a 2 e ~2at + (a *)2 e -2a-t 



4tt 



mf 



The Hubble parameter ( 4.34 ) corresponds to the following scale factor 

1 



a(t) = ciq exp < Hot 



8m 2 



Pa_ -2at 1 Pa* -2a* t 



or 



*2 



(4.36) 
are real 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



substituting the explicit formula for p a and p a * we get 



'(£ -2) £ z e 
a(i) = ao exp <j //q^ — ( — : — o — cos(2^t) 



2 c -2rt 



Am 2 



2m 2 , 



[(r 2 - v 2 ) cos(2i/i) - 2ri/ sin(2i/t)] 



(4.43) 

We see that a late time expansion regime corresponds only to Hq > (compare with pC|]), 
Let us note that a constant part of the Hubble parameter H$ can be incorporated in a 
plane solutions [20]: 

Mt) = A n e a ^ + B n e~ a ^\ (4.44) 

where a n ,H is related with a n as 

a l,H + 3H a n! H = o? n . (4.45) 

Using this approximations authors has been investigate the case of the SFT inspired values 
of parameters £ 2 and c and obtain a cosmic acceleration with a periodic crossing of the 
w = — 1 barrier. In our paper we obtain the similar result for arbitrary values of £ 2 and c. 
Note that for some particular values of these parameters we can obtain one-fold crossing 
of the w = —1 barrier. 
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4.3.2 Two Real Roots Solutions 

The above-mentioned solutions for real roots m n correspond to monotonic behaviour of the 
Hubble parameter. To describe nonmonotonic behaviour let us consider the case of c > 1 
and £ 2 < ^ax- There exist two real roots of ((Of ) mi and vn<i such that that \m\\ < \rri2\- 
The corresponding solution to ( |3,5D is 

tj, = fa + fa, where fa = Ae mit , fa = Be m2 \ (4.46) 

where A and B are constants. Without loss of generality we can put A = 1. 
Using Q4.8] ) we obtain 

ti = ~ 2m| { p ^ e2mit + B2 Pm 2 e 2m2t ) • (4.47) 

Let us analyze a possibility H = 0, which correspond to crossing of the cosmological 
constant barrier for the state parameter w. In Subsection 3.2.2 we have obtained that 
Pmi < and p m2 > (See Figs, [l] and ||). So, for any roots m\ and m2, there exist such 
real B that = at the point t = t\\ 



B= ±J Pmi 9 6 . . (4.48) 



We conclude that solutions (4.46) correspond to cosmological models with the crossing 
of w = — 1 barrier. 

The Hubble parameter and the scale factor are as follows: 

H = - — ^ { Pmi m 2 e 2m ^ + B 2 p m2 m ie 2m2t ) + H , (4.49) 

4mpmim2 ~ 

a = a exp ( - 1 (p mi m^ 2mi * + £ 2 f> m2 m?e 2m2t ) + i? i J • (4-50) 
If mi > and m2 < 0, then at late time 

^~ "2m| Pmie2mii >0 (451) 



and 

fl * " ^•» e2 ™" > °- < 452 > 

Using 

= W(V>i,Va) = 1—5 (Pm 1 m 2 tpl +pm 2 mi^l) +H , (4.53) 

4m^mim2 

we obtain the fourth degree polynomial potential 

V{fa,fa) = 3m 2 (^ 4m 2 mim2 (Pmi m 2^i +Pm 2 rn iV ; 2) ~ ^ - ^ (Pmi^i + Pm 2 ^l) • 

(4.54) 
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We can conclude that all solutions ( [4.46] ) correspond to cosmological models with the 
crossing of w = — 1 barrier. Let us remind in this context that models with a crossing 
of the w = — 1 barrier are a subject of recent studies [|(], 41, 30, 61, 62, 63]. Simplest 
models include two scalar fields (one phantom and one usual field, see [38, 54, |65|] and refs. 
therein). In our case a nonlocality provides a crossing of the w = —1 barrier in spite of 
the presence of only one scalar field. This fact has a simple explanation. The crossing of 
the w = —1 in our case is driven by an equivalence of our nonlocal model to a set of local 
models some of which are ghosts. 

4.4 Cosmological models for N-mode solutions 

Let us construct a cosmological model, with H, defined by ( [4,12| ). To do this we use the 
superpotential method |59| . We consider the Hubble parameter as a function (superpoten- 
tial) of ip n : H(t) = W(ijji,ip2, • • • , ^jv)- To construct potential we use the simplest form of 
the superpotential 



N 



H(t)=WWi,i> 2 ,...,i> N )=J2 W nWn) 



n=l 



Using formula 



H 



^dW n 



71=1 



dlpr, 



— V ( 2ce" 
2ml ^ V 



2«£ + £2 



we obtain, that (4.9) is satisfied if 

dW n = 

In the case of simple roots a n the solutions 

V>„ = A n e a ^ + B n e~ ant 
satisfy the following first order differential equation: 

a 2 n (ipl - AA n B n ) . 



So, 

and 



dWn 
dlpn 



Pc 



A. 



Pa, 



2i7ipa, 



W>2 -4A n £ r , 



■ p an 

2m? 



^fipl - AA n B n - 2A n B n In (a(ip n + sjtfi - 4A n B, 



The corresponding potential V(ipi, . . . , ip n ) is 



(4.55) 
(4.56) 

(4.57) 

(4.58) 

(4.59) 
(4.60) 

(4.61) 



V(Vi,...,V>n) = 3m^ 2 -^|^2 



2a? 



3m r 



H W A +E^(^-4^n). (4.62) 



We see that the potential and does not depend on sign of W . The potential V(ip) is not 
a polynomial, but in the limit of the flat space-time (nip — > oo) we obtain a quadratic 
polynomial. Formula ( [4.62 ) is a straightforward generalization of ( [4.19 ). Note that the 
form of potential is not unique (compare with 
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5. Conclusions 



We have studied linear nonlocal models which violate the NEC. The form of them is 
inspired by the SFT. These models have an infinite number of higher derivative terms and 
are characterized by two positive parameters, £ 2 and c. 

The model with c = 1 is a toy nonlocal model for the dilaton coupling to the gravitation 
field. A distinguished feature of it is the invariance under the shift of the dilaton field to 
a constant. 

For particular cases of the parameters £ 2 and c the corresponding actions describe 
linear approximations to the bosonic [44, |3^[ and nonBPS fermionic j|6| cubic SFT as well 
as to the nonpolynomial SFT [g7|, £||. 

The case £ = corresponds to a linear approximation to the p-adic string (Tf|. Let us 
note that recently a p-adics string inflation model has been considered [23|. 

In the flat case all solutions of the equation of motion are plane waves and are controlled 
by roots of the characteristic equation. Our characteristic equation has complex distinctive 
simple roots. In some particular cases there are single or double roots, which are real or 
pure imaginary. The energy on plane waves is equal to zero except for the cases of couples of 
roots (a, —a). The pressure is a sum of one mode pressures. The pressure for the one plane 
wave corresponding to a real root can be positive or negative depending on parameters of 
the theory. For c < 1 the one mode pressure is positive and for c > 1 it could be negative 
or positive. 

To study cosmological applications we have investigated the behaviour of the models 
in the Friedmann background. We have performed this study within an approximation 
scheme. A simplest approximation is a local field approximation (or a mechanical analo- 
gous model in a terminology of Jl8|]). On an example of the free flat case we have shown 
that in special cases we can use a local two derivatives approximation, but the next deriva- 
tive approximations exhibit artifacts. Followed [22] we have used the Weierstrass product 
representation to study finite mode approximations. As was noted in [22] a straightforward 
application of the Ostrogradski method to these approximations indicates that energies are 
unbounded (an eigenvalue problem for the unbounded hyperbolic Klein-Gordon equation 
on manifolds is solved in |66[| ) and it is expected [^] that an incorporation of non-flat 
metric or nonlinear terms could drastically change the situation. 

A distinguished cosmological property of these models is a crossing of the phantom 
divide |20| , [24]]. But there are also possibilities for other types of behaviour. Namely, 
the toy dilaton model possesses decreasing solutions describing asymptotical flat Universes 
(adding a cosmological constant modifies these solutions to near-de Sitter solutions). It also 
has odd bouncing solutions describing a contracting Universes meanwhile even bouncing 
solutions are forbidden. For special values of parameters corresponding to tachyon SFT 
models there are even bouncing solutions with an accelerated expansion. 

We have shown that for some particular cases there are deformations of the model such 
that exact solutions of the linear problem are inherited by nonlinear non-flat ones. This is 
similar to what was done before for local models [37, [37], ^38|. A stability of exact solutions 
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in local models has been studied in |}7]]. We will study stability of our solutions in the 
future work. 
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A. Calculations of Nonlocal Energy Density and Pressure on Plane Waves 

In this appendix we calculate the energy density and pressure for the following solution 

cj) = Ae ait + Be a2t , (A.l) 
where a\ and a 2 are different roots of (|3.6|). We have 



E{Ae ait + Be a2t ) = E{Ae ait ) + E{Be a2t ) + E cross (Ae ait , Be a2t ) = ABE cross (e ai \ e a2t ), 

(A.2) 

where the functional E cross is defined as follows: 



E, 



cross \ <P1) 92 1 



+ 0l4>2 + C$i$ 2 + 



+ C 



1 r 



e-^) d 2 (e^^) + (e-" d2 <5> 2 ) d 2 (e^^ 

$ 2 = e~ d2 (f) 2 . 



dp 



dp, 



$i^e- 92 h, 



For 4>i = e ait and 4>2 = e" 2 * we have 



E 



cross \€ , 6 



ait p a 2 t) = e («i+«2)t| aia2 _ 1 + ce (-« 2 -«i) + 



+ ce~ a2 - a2 



o 



a 2 e 



ol\u 2 I e 



dp}- 



If ai = — a 2 , then it is easy to show that 



E f P «i* p-arf) — -2n 

■ LJ cross\ c j c j — ^ai 1 



where p a is given by ( pi ). We get 

£(,4e Qli + Be- ait ) = -2AB Pai . 



(A.3) 



(A.4) 



- 28 - 



In the opposite case («i 7^ —0:2) 

E cr0 ss(e ait , e° 2 ') = E aua2 e^ +a ^; (A.5) 

{-a 2 -a 2 ) ( 2 a\-a\ _ 2 a 2 -a 2 



^ ll02 =ai« 2 -l + 5 5 — • ( A - 6 ) 

a 2 — af 

Constants a\ and a 2 are roots of (|3.6| ), therefore, 

= 2 1 2 {«i«2 ((e 2 a| - ce" 2a2 ) - (£ 2 a 2 - ce" 2 " 1 )) + a 2 £ 2 a 2 - a 2 £ 2 a 2 ) = 0. 

a 2 ~~ a l ^ > 

(A.7) 

Note that, the equality E cross (e ait , e" 2 *) = at ai / — «2 also follows from the energy 
conservation low. Let us calculate the pressure P{4>) for the solution <p = Ae ait + Be a2t . 

P (Ae ait + Be a2t ) = P (Ae ait ) + P (Be a2t ) + P cross (Ae ait , Be a2t ) , (A.8) 

where 

Pcross(4>l,<t>2) = ^A erw (01,^)+^2 I ^(01j^)-^Poro«(01»^)-^tJl»BM(01>^2) J ( A - 9 ) 

E kaross = fdfadfo, E Pcross = fafc + C $!$2, (A.10) 

^nzi_ ss = c j [(V^i) d 2 (eT^^) + (e"^) d 2 (e^i)] dp, (A.ll) 

and 

"1 r- 



jf [a (e^Si) 8 (e-^ 2 $ 2 ) + 9 (V^ 2 $ 2 ) d (V pa2 $i)] dp. (A.12) 



The functional E cross = E kcross + E ni2croa3 + E Pcrosa + E n u croaa , proving that £ cross = 0, 
we obtain that 

E k + E nl2 = and E v + E nll = 0. (A. 13) 

Across 1 it><>£cross y cross 1 /it J. cross V / 

Therefore 

Pcross {e ant , e akt ) =0, if (A.14) 
The straightforward calculations give that 

Pcross {e a "\ e~ ant ) =0. (A.15) 

So, P cross (e ant , e akt ) = for all a n and a fc / a n . The pressure is as follows: 

\n=l / n=l 

where 

P(e at ) =;Pae 2at . (A.17) 
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